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Abstract. As a first step towards a successful field theory of Brownian particles in 
interaction, we study exactly the non-interacting case, its combinatorics and its non- 
linear time-reversal symmetry. Even though the particles do not interact, the field 
theory contains an interaction term: the vertex is the hallmark of the original particle 
nature of the gas and it enforces the constraint of a strictly positive density field, as 
opposed to a Gaussian free field. We compute exactly all the n-point density correlation 
functions, determine non-perturbatively the Poissonian nature of the ground state and 
emphasize the futility of any coarse-graining assumption for the derivation of the field 
theory. We finally verify explicitly, on the n-point functions, the fluctuation-dissipation 
theorem implied by the time-reversal symmetry of the action. 

PACS numbers: 05.40.Jc 
Introduction 

Glassy features in the dynamical behavior of colloidal suspensions are among the most 
intriguing aspects of their phenomenology. Their idealization as fluids of interacting 
Brownian particles might provide the key to understand glassiness for continuous, finite 
dimensional systems. A satisfying analytic treatment seems, however, still out of reach. 
The evolution of a system of N Brownian particles can be described either by a system 
of N stochastic ordinary differential equations or, after introducing a density field, by 
a stochastic partial differential equation pQ, [2] which in turn can be recast into a field 
theory. A field-theoretic approach is the ideal framework for the study of conservation 
laws, encoded as symmetries of the action. The appeal of a field theory is then 
twofold: the ease to conceive perturbative series expansions in the (weak) inter-particle 
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potential, and the possibility to obtain non-perturbative results deriving approximate 
self-consistent equations for the correlators. In either case, the symmetries of the action 
should act as a guidance and a constraint for any approximation scheme. 
In order to study the glassy dynamics it is necessary to examine a field-theory in its 
strongly-interacting regime and, to this aim, attempts have been made to re-sum whole 
classes of diagrams, at all orders in the interaction potential [3], [I], [5], [6]. The main 
obstacle towards a successful resummation scheme is that the time-reversal symmetry of 
the appropriate action involves non-linear transformations of the fields, relating n-point 
functions for different n's in an intricate fashion: this precludes the possibility of writing 
easily-solvable self-consistent closure relations for the n-point functions which respect 
the symmetry. In general, a symmetry enforces relations between the field-theoretic 
correlators; in particular, the time-reversal symmetry of the field theory we are about to 
study leads to the correct fluctuation-dissipation relations expected between correlation 
and linear response functions. The fluctuation-dissipation theorem (FDT) is a firmly 
grounded physical requirement for an equilibrated system: solutions violating FDT may 
arise because of terms which drive the system out of equilibrium, explicitly breaking 
time-reversal symmetry, or as a consequence of the spontaneous symmetry breaking 
occurring in glassy phases. However, any solution describing the liquid phase would 
clearly be unsatisfactory if it were incompatible with FDT. 

The aim of this paper is a thorough study of the Brownian gas, which is the th order 
term of any perturbative expansion in the inter-particle potential. Here, the challenge 
posed by the non-linear symmetry is isolated from further complications due to the 
interactions, the origin of the non-linearity being purely entropic. The material is 
organized as follows: in Section [T] we define the Brownian gas and compute exactly the 
space-time dependence of all the n-point correlation functions for the particle density. In 
Section [2] we briefly review how to formulate a path-integral description of the system 
and solve the field theory exactly, re-computing all the n-point correlation functions 
in agreement with the results of Section [TJ Calculations with field-theoretic methods 
are non-trivial since the Martin-Siggia-Rose (MSR) action describes an interacting field 
theory, even though the original system contains non-interacting particles. In Section 
[3] we emphasize the role of the vertex in enforcing the constraint of a strictly positive 
density field as opposed to a Gaussian free field. No coarse-graining procedure and 
no cut-off scale are in principle necessary for a proper definition of the theory. A 
central, non-perturbative result is the proof that the statistics of the ground state of 
the field theory is Poissonian. In Section H] we compute the Jacobian of the non-linear 
field transformation and show explicitly how the n-point functions computed in the 
previous sections satisfy the fluctuation-dissipation theorem associated to the time- 
reversal symmetry of the action. 
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1. The particle formalism 

1.1. Brownian particles 

In natural units (kg = 1 and drag coefficient 7 = 1) the equation of motion for a 
Brownian particle labeled by j moving in M d is: 

x j (t)=rj j (t) (1) 
r]j is a Gaussian white noise with probability distribution: 

P[ Vj ] oce-3W^ 2(i) 
and correlators: 

fo,«(*)%&(0) = 2T6 i:j 5 ab 6(t - t') 
where a and b label the vector components. The solution of equation ([I]) is: 

x j (t)=x 3 (0)+ [ dt' Vj (t') (2) 
Jo 

Note that for any j, Aj(t) = Xj(t) — Xj(0) = J dt'rjj(t') is a Gaussian random variable. 

Soon we will be interested in the value of {(Ax a (ti) — Ax b (t2)) 2 ) v , which can be easily 
computed using: 

(Axi^Ax^ih)^ = [ dt' I dt" {m<Jf)rijjiSf))v 

Jo Jo 

= f 1 dt' f 2 dt" 2T5 lj 5 ab 5(t' - t") = 2T5 ij 5 ab mm(t u t 2 ) (3) 
Jo Jo 

The result in ([3]) is non-vanishing only if Ax's relative to the same particle are considered 
and finally: 

((Ax a (ti) - Ax b (t 2 )) 2 ) v = 2T5 ab \t 1 - t 2 \ (4) 

1.2. The Brownian gas 

The Brownian gas consists of N Brownian particles ([1]) confined in a (i-dimensional box 
of volume V with periodic boundary conditions. As usual, the thermodynamic limit is 
defined as: 

N 

V — » 00 ; A^^oo ; — po = const (5) 

The free parameters describing the gas are its average density po an d its temperature 
T, or equivalently the diffusion coefficient related to the temperature by the Einstein 
relation: D = In natural units 7 = 1 and kg = 1, so that D = T. The 

natural length-scale in the system is the average separation between particles, defined as 
A = p 1 ^ d . Also, the well-known relation for diffusive processes (Ax 2 ) = 2D At = 2TAt 
fixes the natural time-scale of the system to be r = Once distances and time 
intervals are expressed in these natural units, the Brownian gas becomes a parameter- 
free universal system. 
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1.3. Density correlations 

The number density field g(x, t) is defined as: 

N N 

g(x,t) = Y,S(x-x 3 (t)) or: g k {t) = -^e lk ^ (6) 

3=1 3=1 

We now compute the correlation of densities measured at n different points in space and 
time. The average over initial conditions and realizations of the noise is fully determin- 
istic and represents, at least at equilibrium, the expected outcome of a measurement. 
The average over the initial positions (. . )i.p. and the realizations of the noise (. . .) v are 



defined in Appendix A 



The general structure of the following computations is: 

- The position of each particle at time t depends both on the initial position and 
on the realization of the associated noise. Using (J2j) the two contributions can be 
separated and the averages over initial positions and noise realizations factorize. 

- The average over initial conditions is computed using formula (1A.1I) . 

- The average over the noise is performed exploiting the standard relation (e~ 1 ^)^ = 
e -5<V ,2 > 1 /> ; valid for any Gaussian random variable ip with probability distribution 
P [■?/>] oc e - ^ 2 , and applying ([3]) and (jlj). 

The final result for a generic connected n-point function is surprisingly compact and is 
reported in formula (jTUI) . 

1.3.1. 2-point functions It is a simple exercise to verify that the 1-point function yields 
Po, as expected from ([5]). Applying the steps described above, for the 2-point function 
we obtain: 

(gix^tM^h))^. = £ (eM^)^ e -r^-^ + 

k 

The connected correlation function is: 



V 2 \^ J V 2 



(g(x 1 ,t 1 )g(x 2 ,t 2 )) c = (g(x 1 ,t 1 )g(x 2 ,t 2 )) ri ,i.p. - {g{.x l ,t 1 )) v j. P Xg(x 2 ,t 2 )) ri j.p. 

k 

In the continuum limit another factor of V appears in the numerator of the first term 



k 



k k 

which is then finite in the thermodynamic limit whereas the second term vanishes, being 
suppressed by a factor y. In the thermodynamic limit the 2-point function is thus: 

d d k 



(^i,ti)^(x 2 ,t 2 )) c ^p / e ^H^) e -Tk^t 2 \ 



1 f (x 1 - x 2 



(47rT|t 1 -t 2 |) rf / 2 v \ AT^-t 



\2 
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1.3.2. 3-point functions The computation of the 3-point function gives some further 
useful insight. The generalization to an n-point function is then straightforward. 

(g(x 1 , t 1 )g(x 2 , t 2 )g{x 3 , t 3 )) V:LP . 
= e ikl - Xl e ik2 - X2 e ik3 - X3 ^5 kl+k2+k3 N(e-' lkl - A ^ 



V 3 



+ e tk2 ' X2 e ik3 - X3 ^S k2+k3 N(N - i)( e -^2-A(t 2 ) e -*fc 3 -A(*3)^ 
+ (1 <-> 2) + (1 <-> 3)) + -^-N(N -1)(N- 2) 



V 3 

Each term in the calculation above is associated to a connected or a disconnected 
correlation; connected pieces represent single-particle self-correlations to which a formal 
"conservation of momentum" rule is associated by formula (lA.lj) . Again, the subtraction 



of disconnected pieces removes, in the thermodynamic limit, all the terms except for the 
first one, which is the only fully connected term. 

1.3.3. n-point functions A connected n-point correlation functions can then be written 
in the thermodynamic limit as: 

{g{x 1 ,t 1 ) . . . g(x n ,t n )) c (8) 
= e ikl - Xl ...e ik "- x "— 5 kl+ ... +kn (g{ki,ti)...g{k n ,t n )) v , L p. 

where (g(k\, h) . . . g(k n , t n )) v = ( e -' fc i- A (ti) . . . e ~ lkn ' A( - tn ^) v . After computing the average 
over the noise, the general n-point correlation function is easily obtained: 

n 

(gih,^) . . . g(k n ,t n )) v = exp (T^kj ■ kj\U - tM (9) 

i<j 

Note that S kl+ ,,, +kn in ([8]) allows to perform easily one of the sums over k and the other 
n— 1 summations give, in the continuum approximation, a factor V n ~ x as in (J7J). Hence 
the whole expression is proportional to po and it is finite in the thermodynamic limit: 

(g(x 1 ,t 1 ) . . . g(x n ,t n )) c (10) 

= [ d ki r d fc n _i jk v (xi-x n ) ik n . 1 -(x„. 1 -x n ) TY:7 < j k i- k i\ t i- t j\ 
H0 J (2ir) d "'J (2vr) d 

with the constraint k n = —(k\ + . . . + fc n _i) in the last exponential. The space and 

time-translation invariance of the result are manifest. 

As we have seen, the only finite contributions in the thermodynamic limit come from 
single-particle self-correlations. The particles' positions Xj(t) obey the Langevin equa- 
tion ([1]) and the solutions are Gaussian random variables for which all the connected 
moments (x(ti) . . . x(t n )) c vanish if n > 2. It could be somewhat surprising to find 
non-vanishing connected moments of the density for every n. The key point is that the 
density field is a function of the particle positions and a generic function of Gaussian 
random variables does not have to be, and in this case is not, Gaussianly distributed. 
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2. The field-theoretic formalism 

We now briefly review how to introduce a single differential equation with all the 
information content of the system of Langevin equations describing the Brownian gas, 
and how to rephrase it as a field theory. 

2.1. Dean's equation 

The Brownian dynamics of a system of particles can be encoded in a stochastic ordinary 
differential equation for each particle; since we consider the non-interacting case only, 
the equations are naturally decoupled: 

ij(t)=Vj(t) ; j = l,...,N 

Defining the density field as in ([6]), the system of equations can be formally translated 
into Dean's equation ( TTTT) . which is the corresponding partial differential equation 
governing the evolution of the density field: 

d t Q = V-(^) + TV 2 g (11) 

In the process, the N independent noises rjj(t) have been turned into a single vector 
field £ (x, t) with self-correlation: 

(£ (z, t)£ b (x', t')) = 2T5 ab 5(t - t')5(x - x') 

Since in ( TTTT) the noise multiplies the density field, the differential equation becomes 
a meaningful string of symbols only once a proper interpretation for the noise term is 
given [?]. We adopt the Ito prescription, which implies a trivial Jacobian [I] for the 
path- integral in ffT3l) . 

It should be stressed that, due to the singular nature of the density field, the natu- 
ral setting of Dean's equation is a distribution space. In such a setting the equation has 
been derived exactly, without any coarse-graining assumptions or cut-off scales which 
are by no means necessary to obtain the differential equation [2]. 

2.2. The MSR action and the Feynman rules 

2.2.1. The MSR procedure Through a formal procedure first devised by Janssen [8] 
following Martin, Siggia and Rose (MSR) [9], a stochastic differential equation can be 
recast into a field theory defined by a functional integral. The basic idea is to impose 
the original equation of motion for the density field p by introducing a new field i(j) as 
Lagrange multiplier, whereas the noise is integrated out exploiting its Gaussian measure. 

After introducing for convenience a shifted density field with a vanishing average 0, 

g(x , t) = p + p(x, t) s.t. (g(x,t)) = p , (p(x,t)) = (12) 

| The spatial average is always well defined, being one of the parameters defining the thermodynamic 
limit. At equilibrium we expect it to be equivalent to a temporal average at a fixed position and also 
to an average over the noise histories at fixed position and time. 
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the analytic expression of the MSR generating functional is: 

Z = J VcJ)Vp e s (13) 

with action: 

S = J d d x dt {i4>[d tP - TV 2 p] + Tp [Vi4>} 2 + Tp[Vi0] 2 } (14) 

The time integration in ( Tl4l) extends from — oo to +oo or, equivalently, the initial 
conditions are shifted to the infinite past: due to the dissipative dynamics, the ground- 
state of the theory describes then the equilibrated system, forgetful of any initial 
condition. We do not specify instead the functional domain over which the path integrals 
in fflBl are computed and proceed formally: the functional integrals are regarded 
as formal series of Gaussian integrals computed over unconstrained (i.e. positive and 
negative) fields whereas, in principle, the integral over p should be restricted to sums 
of <S's only, as the form of the density field in (J6]) would suggest. The exact match of 
the results we will obtain with the ones calculated in the particle formalism justifies 
a posteriori such an approach, emphasizing the role of the auxiliary field icf) and of 
the vertex in imposing the equation of motion (fTTj) and enforcing the positivity of the 
density field. 



2. 2. 2. Feynman rules For the following analysis it will be convenient to work in a rep- 
resentation involving times t and momenta k. Introducing the time-ordering convention 
to write earlier fields to the right of later fields, the appropriate Feynman rules are then: 

(pp)(k,t)=p e- Tk2 \ t \ 

(i}p)(k,t) = 6(-t)e +Tk2t 

{fd})(k,t) = 6(t)e- Tk2t 

(i4>i<f>)(k,t) = 

-2Tfci • k 2 



The dotted edges on the vertex indicate the action of spatial derivatives. All the mo- 
menta are considered positive when incoming in the vertex. 

The propagator (pp) is the density-density correlator and we will refer to the propa- 
gators (icf>p) and (picfi) as the advanced and retarded response functions respectively. It 
is worth emphasizing that such terminology would be appropriate only in the case of 
Langevin equations with additive noise, for which the occurrence of icf) in a field-theoretic 
expectation value can be proven to represent the linear response to the variation of a 
chemical potential p coupled to the density field. For the MSR field theory associated 
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to Dean's equation, instead, response functions are characterized not by i<p but by the 
composite operator V ■ (gVi<f)) as we will show in Section H~T1 Yet, for simplicity, in the 
following we will refer to the pseudo-response functions (i(f>p) and (pi(f>) as "response 
functions" tout court. 



2.2.3. Diagrammatics Although corresponding to a non-interacting particle system, 
the action ( Tl4"l) has a vertex and the field theory is non-Gaussian. The study of the gas 
is nevertheless easier than the interacting case since all loop diagrams vanish because of 
causality, as the following example, easily generalizable to arbitrary n-point functions, 
illustrates: 




Labeling t 1? t 2 and £3 the vertices of the triangle, the ^-functions in the propagators im- 
pose: t\ < t2, ti < t% and t 3 < t\. Such constraints can be satisfied only by ti = t 2 = £3, 
with a vanishing volume in phase space. Hence, only tree diagrams yield non-vanishing 
amplitudes and such a simplification which will allow to compute exactly all the corre- 
lation functions. 



The structure of the vertex implies that among the edges of a tree there is one and 
only one correlation-like propagator, all the other propagators being response-like. The 
identity below allows to express a correlation-like propagator as a sum of two response- 
like propagators: 



(pp)(k,t) = p e- Tfc2|t| 



Po 



Po 



6(t)e 



-Tk 2 t 



6(-t)e 



TkH 



(pi<f)){k,t) + (i(f>p)(k,t) 



(15) 



Graphically: 



Po{ + )• 

each diagram can be decomposed into the sum of two simpler di- 



Using relation ([151) 

agrams in which all the propagators are response-like: causality, encoded through the 
6{t) functions in the Feynman rules, endows the response-like propagators with a natural 
orientation so that the concept of "earlier vertex" is well-defined. The above decompo- 
sition leads also to the appearance of vertices with three dashed legs whose properties 
are discussed in Section 12.41 



After applying the decomposition ( TTBI) the diagrams contributing to a correlation func- 
tion can be divided into two subsets: 



- c-diagrams: all the external lines are solid lines or, equivalently, one of the vertices 
has three dashed legs. 
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- r-diagrams: one of the external lines is a dashed line or, equivalently, all the vertices 
are of the kind depicted in the Feynman rules of Section 12.2.21 

Indicating by C and R the sum of all the c-diagrams and r-diagrams respectively, we 
have: 

- Response functions = R 

- Correlation functions = C+R 

In the next section we will develop a computational scheme for calculating all the 
diagrams of the theory and, after summing together the appropriate contributions, 
we will give explicit analytic expressions for the generic n-point density correlations. 
In order to simplify the notation, the temperature will be set to one in the following 
(T = 1) an it will be easily reintroduced in the final result (135]) . 



2.3. Computation of all the r-diagrams 

r-diagrams can be immediately mapped to binary rooted trees by considering the prop- 
agator with the external dashed line as the root of the tree. Before confronting the 
general case it is instructive to compute explicitly the first few r-diagrams to better 
understand and appreciate the combinatorial approach developed in section 12.3.41 The 
hasty reader might want to continue there directly. 



Conventions: 

ti < ti < . . . < t n _i 

t n 

m ab...yz '■= m l n (ta, tb, ■ ■ ■ , ty, t z ) 

2.3.1. 3-point Applying the Feynman rules and constraining the integration domain 
according to causality: 

1 



2 3 

R 3 (l, 2, 3) = / dt {-2k 2 ■ fc 3 ) e - fc ?(*-*O e -*2(fe-*) e -k§(ts-i) (16) 
Jti 

where conservation of momentum: — k\ + k\ + h\ = —2k 2 ■ h% has been used. 

As explained below, the result can be depicted as: 

1 

12 A 1 



R 3 = /\ - /\ (17) 
^ 2 3 2 3 v ; 

In the following we will be careful with the terminology, calling "diagrams" the usual 

Feynman diagrams [10], mapped to integrals by the Feynman rules, and "graphs" the 

objects of the kind drawn in ( flTI) . which are a graphical representation of the result of 
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such integrals. 

The crucial property of such a graphical representation is that the group of permutation 
of labels under which the exponentials in an analytic term are invariant is isomorphic to 
the group of topologically equivalent relabelings of the associated graph. As we will see 
starting from the computation of the 4-point r-function, the graphical approach greatly 
eases the bookkeeping. 

The rules for drawing a graph given the result of the integrations are very simple: 

- To each term of the form e~ k ^ a ~ b " > associate two vertices labeled 'a' and l V linked by 
a line. The edge should be labeled by yet we will avoid doing so since the labels 
for the edges of a graph can be easily guessed considering the momenta labelling 
the original Feynman diagram and the contractions to obtain the graph. 

- Identify all the vertices with the same label 

Note that no loops are generated as a consequence of the tree-like structure of the 
original Feynman diagrams and tadpoles are naturally removed since a = b implies 

Recovering an analytic expression given a labeled graphs is an obvious procedure. 
Our graphical convention is to orient both, diagrams and graphs, in such a way that 
the arrow of time points towards the bottom of the page. Then, the label m a i, yz as- 
sociated to the vertex of a graph is the time labeling the earliest of the underlying leaves. 



2.3.2. 4-point The 4-point r-function can be easily computed from the result of the 
3-point r-function: 



r 4 (l, 2, 3, 4) = / dt (+2k 2 ■ k^e-^-^e-^^R^A, 3, 4) 
Jti 

-kj(m234-ti)-k%(t2-m234)-k\(m 3 4-m 2 34)-k%(t' i ~m' i 4)-kl(t4-rn 3 4) 



_ e -ftnm234-tlj-fc5(t2-m234j-KAl m 34-m234j-K3(,ta-m 3 4 / )-fc|(t4-m34 - ) 
_ g~ fe 2 (*2 _ *i )~k\ ( m 34 — ti ) — fe| (*3 —mai)—k\ {U -m.34) 

-k\ (m-234 ~t\ )-&2 (*2 -m 2 34 )— fc§ (*3 -m 2 34 )-k\ (*4 -m-234 ) 



( -\- k 2 ~\~ k^ -\- k^) 

_\ 2 ^ 2 • ^ A e -fc|(t2-tl)-fci(t3-tl)-fc!(t4-tl) 

( k-y -\- ~\~ k^ -\- k^) 

where &u = k\ + k 2 = —(k 3 + k A ). In the last two terms of the result, rational 
functions of the external momenta make their first appearance. The exponentials and 
the denominators in the last two lines are independent of the permutation of labels 
{2, 3, 4} and since the interest is eventually in the sum over all the topologically distinct 
permutations, we add together the numerators of the rational expressions; due care 
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should be taken since k A does depend on the particular permutation and the k A s relative 
to different permutations need to be distinguished: 

2(k 2 ■ k A + k 3 - k A > + k 4 ■ k A i>) 

= 2{k 2 ■ (kt + k 2 ) + k 3 ■ (hi + h) + k 4 ■ {h + k A )) (19) 
= 2{k\ + k 2 3 + kj + h(k 2 + k 3 + h)) = 2(-k\ + k 2 2 + k 2 3 + k\) 

Simplifying the denominators, both rational expressions disappear, leaving a factor of 2. 
The final result has then the following representation, the capital 'R' in i? 4 representing 
the sum over all the possible labelings of all the r-diagrams: 




TI 



Y2ti ^ s a sum over a ^ topologically inequivalent relabelings of the graphs. Note 
that the property of two labelings of being TI (topologically inequivalent) or TE 
(topologically equivalent) depends on the particular graph. The first two graphs in ( f20l ) 
need to be summed over the 3 TI permutations of external labels {2,3,4} ('1' being 
always fixed since any r-diagram vanishes unless the external dashed line is associated 
to the earliest time). The last two graphs do not need any further treatment or, 
equivalently, the 3 possible labelings are topologically equivalent and they have been 
summed over already to replace the rational functions with the integer '2'. 



2.3.3. 5-point In this example one last complication, given by the existence of several 
topologies of binary trees with the same number of external legs, needs to be faced. In 
order to remove the rational functions it is required to sum not only over the permu- 
tations of labels, as before, but also over the possible topologies of trees. Note that in 
this context diagrams with the same topology but whose propagators have different ori- 
entations are considered topologically distinct and they are mapped to different binary 
rooted trees, namely 5a and 5b below. 



Topology a: 

1 




r "12345 

r 5a (l, 2, 3, 4, 5) = / dt (+2k 2 ■ k B )e- k ^-^e- k ^-^r 4 (B, 3, 4, 5) = 
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A 2 2 



i 

[2 

3 4 5 



,3 
3 4 5 



+ A 2 2 /l\ + A 3 



2 3 4 



2 34 5 



where: 



Ai 



A^ = A 



2k 2 ■ k B 



A, 



2k 3 ■ k A 



Topology b: 



k^ 4- k 2 ~\~ k 3 4- k^ 
2k 2 ■ k B 

2 7 2 7 2 7 2 7 2 7 



°y/\ A 



2 3 4 5 



~kg ~\- fcg -|- -|- fcg 



(22) 



/■m.2345 

r 5b (5, 4, 3, 2,1)= / rft(-2fc A • fc^e"^"'^^, 2, 3)i? 3 (A 4, 5) 



1 

v 2 
2/ \4 



2/ \4 
2 34 5 2 34 5 

1 

2 



1 

23 A + Bj 2 ¥\ 

4 5 4 5 



(23) 



2 3 



,1 

\ 2 
2 3 



B 2 4 V\ + £ 3 



2 3 4 5 



B 3 



2 3 4 5 



where: 

» _ -2fc,4 • fee _ -2fc A ■ k c 

~ -k 2 4- A- 2 4- P 4- P -P 4- P 4- P 4- P 

»vj t "'2 ' 3 ^ A 1 ' 4 ^ 5 ^ C 

B 3 = ~ 2kA ' kc (24) 

k\ 4- 4- fcg 4- 4- k§ 

In the following, with 5 a n th and 5 b n th we will refer to the n t/l graph appearing in (|2~TI) 
and (123]) respectively. 

A subtlety: 

The group Q^ a of TI relabelings for Feynman diagram 5a contains 12 elements and in- 
duces 12 relabelings for the graphs in (j2ip . For each graph in (121]) we can factorize £?5 a as 
a product of TI and TE relabelings, the factorization depending on the particular graph 
considered. The factorization is possible since the graphs generated by contracting edges 
in diagram 5a have a higher symmetry than the diagram itself. Unfortunately this is not 
always the case; in fact, Feynman diagram 5b has only 3 TI relabelings. However some 
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of the graphs in (1231) . namely 5 63 rd, 5& 4«,, 5 65 rd and 5 ft 6 th, admit 6 TI relabelings since 
the edge contractions have lowered the symmetry of 5b. Nevertheless, the troublesome 
graphs appear in pairs ({5 6 3 rd, 5 bj5 th} and {5 6)4 rd, 5 6 6 th}) so that C/5& induces 3 relabelings 
on each of them for a total of 6 relabelings associated to each topology. The 3 elements 
in Gab are thus enough to induce all the TI relabelings for any graph in ( 1231) . In general, 
whenever the contractions reduce the symmetry factor of a diagram by a factor n, n 
copies of the same graph appear, so that the TI labelings of the original diagram are 
sufficient to induce all the possible TI labelings for the topology of the graph. This 
is a general phenomenon deriving from the fact that the graphs are decorations of the 
underlying tree structure determined by the Feynman diagrams. 

The results for diagrams 5a and 5b (( 12T1) and ( 1231) respectively) contain several 
graphs with the same topology: examples of such terms are {5 a3 rd,5 b3 rd,5 b5 th}, 
{5 a 4 th, 5 b4 th, 5 b6 th}, {5 a 7 th,5 b 7 th} and {5 o 8 «», 5 6>8 th}. A labeled graph fully determines 
the exponentials of the analytic expressions, which are invariant under the TE relabel- 
ings of the graph. Instead, the rational functions in ([221) and (T24|) do depend on the 
labeling. As for the 4-point function in ( fl9l) . by summing over all the TE relabelings for 
each graph we now replace the rational functions in (12Tj) and (123"]) with integer numbers. 
The computations are straightforward and can be summarized as follows: for the topol- 
ogy of the groups {5 a3 rd, 5 6 3 rd, 5 6 5 th} and {5 a4 «i, 5 64 t/>, 5 6 6 th} the Feynman diagram 5a 
induces 12 labelings and 5b induces 6. The sum over the 18 labelings yields a factor 
of 2. For the topology of the groups {5 a jth,5 b7 th} and {5 ai8 th, 5 6)8 tfc} the sum over the 
12+3 induced labelings yields a factor of 6. Finally, the result can be expressed in terms 
of a linear combination of graphs with integer coefficients, summed over TI relabelings 
only: 




2.3.4- n-point The previous sections clearly show that the amount of combinatorics 
required to compute higher order r-functions rapidly becomes prohibitive and we need to 
automate the computation of n-point diagrams. This can be achieved in two steps: first, 
starting from a Feynman diagram, generate all the topologies of the relevant graphs, and 
second, associate to each graph the appropriate combinatorial integer factor. Indeed, 
given an n-point Feynman r-diagram, X, there is no need to apply the Feynman rules, 
perform explicitly the computation and then draw the corresponding graphs. Instead, 
the graphs representing the result are readily obtained by the following procedure: 

- Consider all the graphs drawn after contracting m (m = 0, . . . , n — 2) among the 
contractible edges of X, where a "contractible edge" is defined to be either an 
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internal edge or the external edge labeled by the earliest time (which in the r- 
graphs is always drawn as the uppermost edge, with the appearance a stalk). 

- Multiply each graph by a factor (— l) m . 

- To each vertex with I dangling lines associate a factor TZi = (I — 1)!. The integer 
factor associated to the graph is then the product of all the coefficients associated 
to its vertices. When considering an explicit labeling, the label associated to a 
branching is the smallest (earliest) label of its descendant leaves. 

- Sum each graph over its topologically inequivalent (TI) labelings. 

The first two rules follow naturally from an iterative application of the fundamental 
theorem of calculus: f^dtf(t) = F(b) — F(a). In fact, an integral over time is associ- 
ated to each vertex of a Feynman diagram and the evaluation of the primitive at the 
integration boundaries determines the sign and the survival or the contraction of the 
edge in a graph. The third prescription depends on the specific form of the Feynman 
rules and its proof is reported in Appendix B.l All the above prescriptions should be 
familiar after working out explicitly the analytic computation of the 3, 4, and 5-point 
functions of the previous sections. 



The next section aims at motivating one further prescription and at showing explicitly 
the remaining steps for the calculation of density correlation functions. Some readers 
might want to go directly to Section 12.4.41 

2.4- Computation of all the c-diagrams 

In 12.2.31 we have defined the c-diagrams and noticed the appearance of a new vertex 
with three dashed legs which was not explicitly present in the action (fl4l) . We recall 
that the new vertex is a consequence of the application of the decomposition (Tl5|) . After 
summing over the three possible ways of orienting the original r- vertices, or, equivalently, 
of distributing derivatives marked as dots on the edges, it is convenient to introduce such 
vertex, called c-vertex, as a new effective Feynman rule: 

.••'"*••. = *\ + , -'\ + / 

with analytic expression: 

- 2k x ■ k 2 - 2h ■ k 3 ~ 2h ■ h = 2k* + 2k\ + 2k x ■ k 2 = k\ + k\ + k\ (26) 

It is easy to prove that the Feynman diagrams generated with the original Feynman 
rules of section (I2.2.2P and the decomposition (fl5|) are the same as the ones generated 
with a new set of Feynman rules consisting of an r- vertex, a c-vertex and response-like 
propagators only. Diagrams generated only with r-vertices are r-diagrams; the ones in 
which one c-vertex is present are c-diagrams. 

The application of the Feynman rules is straightforward and the computations of the 
integrals very similar to the ones for the r-functions. We just report some results, 
emphasizing the differences with the previous calculations. 
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2.4.I. 3-point 

12 3 

/mi23 
dt (k\ + k\ + jfe2) e -*?(^-«) e -*2 a (te-«) e -*i(t3-«) 
-00 

_ e -fcj ! (ti-mi23)-fc|(t2-nii23)-fc3(t3-nii23) (27) 

Note the different integration domain in comparison with (fl~6]) . The graphical 
representation of (|27l) is: 

C 3 (l,2,3) = /t\ 
v 7 12 3 

The ♦ which marks the c-vertex is just a bookkeeping symbol whose importance will be 
clear soon. 

We are eventually interested in the computation of the 3-point connected density cor- 
relation function which is the sum of the 3-point r-function and c-function. Indicating 
as usual with capital letters the sum over topologically inequivalent relabelings (in this 
case just one), the graphical result is: 

.1 

.'2 



3 



{Q{Kt x )Q{h,h)Q{h,h)) c = R 3 + C 3 = + 

where we have crossed the edges to be removed when applying the contraction rule for 
edges with the same label, as described in section 12.3.11 Finally, the last two graphs 
cancel each other as can be explicitly noted comparing (|T5|) and f[2"T|) . the ♦ being 
irrelevant in this context. In the result, the vertices of the ladder graph are pinpointed 
for clarity by black dots. 

2.4.2. 4-P°^ n t Other n-point functions with n > 3 can be generated recursively from 
one c-vertex and appropriate combinations of r-functions. For the 4-point c-function we 
have: 

/■■\A 

/ I \ 
1 2 ."■ 



/m.1234 
dt (kl + k\ + A^e-^-'V^-^A 3, 4) 
-00 

— e - fc l(*l- m 1234)-fci(t2-mi234)-fcj4(m34-nil234)-fc3(t3-m34)-fc2(i4-'m34) ^£ 



{k\ ~\~ k 2 ~\~ k^ -\- k^ 



4 

or, graphically: 



4^ - c 



c 4 (l,2,3,4) = 23 A - C /A\ (29) 
v ' ' ' ' 45 2345 v/ 
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where: 



C 



(30) 



(kf + k 2 > + kj + kj) 

The Feynman diagram admits 6 TI permutations: differently from the case of the 4- 
point r-function, the allowed permutations for the 4-point c-function involve all the 
labels {1,2,3,4}. The exponential in the second term of (128]) is invariant under such 
relabelings (or, equivalently, the relabelings are TE for the second graph in fl29|) ). 
Since the denominator of C is also invariant, we can just add together the numerators 
corresponding to different labelings. Rewriting (kf + k\ + k\) as (2kf + 2k\ + 2ki ■ k 2 ) 
and summing over the 6 relabelings: 

(2k\ + 2k\ + 2k 3 • k A ) + {2k\ + 2k\ + 2k x ■ k 2 ) + {2k\ + 2k\ + 2k x ■ k 3 ) + 
+ (2kl + 2k\ + 2k x ■ fc 4 ) + (2k\ + 2k\ + 2A; 2 ■ fc 3 ) + (2^ + 2/c| + 2k 2 ■ fc 4 ) = 

— 5(fc^ -|- "I - &3 "I - ^4) 

After the sum, the rational function is replaced by a factor of 5. It is clear that the 
new vertex ♦ generates combinatorial factors different from what we calculated for the 
r-functions and this justifies a new graphical symbol for it. 
The final result for the 4-point c-function is: 



c, 



E 

TI 



(31) 



In order to obtain the 4-point connected density correlation function we add together 
dm and (ED: 



(g(kx, ti)g(k 2 , h)g(k 3 , h)g(k A , t A )) c 



tl 

2 

+ 2 ) ! ( - 5 All = U 




The reader is invited to perform the simple sum over the TI relabelings and verify 
explicitly that, the result for the density correlator is again a ladder graph. 



2.4-3. 5-point The general mechanism should be clear by now. As before, for the 
5-point functions two topologies are present. We only report the final result, after 
summing over the two topologies and all the TE relabelings for each graph: 



C — C*5a + C; 



E 

TI 



+ 



5^-2 



+ 26 



(32) 



Note that the factor of 2 associated to the fourth term of the r.h.s. is a (3 — 1)! from 
the sub-tree with 3 dangling lines. Moreover, an interesting pattern 1, 5, 26, . . . appears 
for the integers associated to a c-vertex with I dangling lines. In Appendix B.2 it will 
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be identified with the sequence of generalized Stirling numbers (Sloane A001705, [TT]). 
Verifying that the sum of fl25l) and fl32l) yields the 5-vertex ladder graph is a lengthy but 



rewarding exercise. In Appendix C.l it is proven for general n-point functions that such 
ladder graphs are the graphical representation of the very same analytic result obtained 
in ([9]) with the particle formalism. 



2.4-4- n-point The general recipe of section 12.3.41 for generating graphs with their 
integer coefficient can now be complemented with a new prescription: 

- To each c- vertex with I dangling lines associate a factor Q = S1-2, where Si is the 
I th generalized Stirling number (see Appendix B.2 for a proof). The c- vertex is 



labeled by '1', being always the earliest vertex in a Feynman diagram. 



2.5. Graph cancellations 

Up to this point we have computed all the r-diagrams and all the c-diagrams writing 
them as sums over topologically inequivalent labeling of graphs. At least in the simplest 
examples, after combining them to compute an n-point density correlation function, we 
have shown that the result is extremely simple, namely a ladder graph, i.e. a single 
time-ordered labeling of a single graph, as displayed in ( [33il . The equivalence of the 
analytic expression corresponding to the ladder graphs and the result for the n-point 
density correlation function computed with the particle formalism in ([9]) is explicitly 



shown in Appendix C.l 



1 

u 





(33) 



Indeed, for every n-point function the same cancellations seen in the examples above 
occur, causing the sum of all the graphs generated from the r-functions and the c- 
functions to vanish, with the only exception of ladder graphs. Let us partition the 
graphs generated by the graphical procedure into three classes: r-graphs with a stalk, 
r-graphs without a stalk and c-graphs. In the figure below we report the example of the 
5-point correlation function: 




+ 6 



+ 2 



(34) 



- 5 



- 2 



26 



/As 
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It can be proven that r-graphs without a stalk completely cancel the c-graphs 
(Appendix C.3) and that r-graphs with a stalk cancel one another with the exception 
of the ladder graph ( Appendix C.2 ). These facts lead us to the following grand result, 
sealing the identity between the particle (Q and field-theoretic approaches: 

■■■g(k n ,t n )) c = tn = p S(k 1 + ... + k n )exp[T^2ki- kjlti-tjl) (35) 

i<j 

where the last term has been complemented with a factor p (2ir) d S(ki + . . . + k n ), po 



coming from the decomposition ( 1151) and the 5 originating from the global conservation 
of momentum associated to the Feynman diagrams. 



3. The Poissonian distribution and a remark on coarse-graining 

After Fourier-transforming (I35I) . the analytic expression for an n-point density correla- 
tion function in direct space is identical to ( flOl) : it follows then that the equal-time limit 
of an n-point function is poS(xi — x n ) . . . 5(x n _i — x n ). At equal times the density field 
is thus delta-correlated in space, as expected for a non-relativistic theory which is both 
causal and local. 



Let us consider the probability distribution function (PDF) for the equilibrium density 
measured at a single point in space-time: the equal-time limit of an n-point function 
is the n th moment of such a PDF. Since a probability distribution is uniquely deter- 
mined by its moments and we know from (1351) that all the cumulants are equal to po, 
it follows immediately that the PDF of the local density is a Poissonian distribution 
with average equal to po- Whereas this statement is obvious when considered from a 
particle viewpoint, it is highly non-trivial from a field-theoretic perspective, involving 
the resummation of an infinite class of Feynman diagrams. 

For an infinite system, A = p is the only natural length-scale. In principle we 
may want to introduce by hand a second length-scale /, which could be a standard unit 
of length (e.g. 1 cm) or being imposed by the spatial resolution of an experiment. Then, 
Po can be considered as a free parameter. In the large po limit the Poisson distribution 
is well approximated by a Gaussian with average po and standard deviation a = ^/po■ 
More precisely, when measured at macroscopic space-time intervals, the 1-point function 
is po ^ 1, the 2-point functions are 0(1) and all higher n-point function are suppressed, 
being Oil/ p^~ 2 ). Hence, the effective behavior of the system can be approximated by 
a free field theory. 

Yet, for po < 1, the Poisson distribution is markedly different from a Gaussian^]: the 

§ No inter-particle potential would ever give rise to a genuine free field theory: in that case, in fact, 
the PDF for the density would be a Gaussian and the density field would not be positive definite. 
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hypothetical experiment is probing the system down to length-scales at which the gran- 
ularity of the gas is manifest. The vertex appearing in the field theory, which faithfully 
reproduces the Poissonian statistics, is then the vestige of the particle nature of the 
gas. Its relevance is crucial at low densities where the free particle system is described 
by a strongly interacting field theory. Note that the length-scale A merely identifies 
the cross-over between "free" and interacting theory but by no means it represents a 
hard cutoff scale for the system. Again, any coarse-graining procedure would be purely 
arbitrary, the field theory being well-defined at any length-scale. 



4. Time-reversal symmetry and fluctuation-dissipation theorem 

This last section contains a brief review of the time-reversal symmetry of the action and 
the fluctuation-dissipation relation it implies. Also, we will finally focus on the mecha- 
nism through which the symmetry is enforced at the level of correlation functions. 

The path integral ( 1131) and the action (114p define the field theory for the Brownian 
gas. The corresponding equation of motion for the physical field g is Dean's equation 
( ITT]) and it can be interpreted as a continuity equation for the density: 

d t g — V • j with current: j = —(,y/g — TV g 



4-1. Response functions 

An external force f(x) acting on the density field shifts the current: 

j(x) -> j(x) +p(x)f(x) = j(x) - p(x)V/i(x) 

the last equality being meaningful only if / is exact and a scalar potential /i can be 
defined: f(x) = —Vji{x). Dean's equation can then be written as: 

d t g = V- {iyfo) + TV 2 g + V • (gVfj) 



,6F 
5g 

where the free energy functional F is defined as: 



d d x 



Tg(x, t) ( In - 1 ) + n(x, t)g(x, t) 



Po 



The physical interpretation of /i as a chemical potential field linearly coupled to the 
density is now fully transparent. Notably, however, in the field theory fi appears in a 
term of the MSR action which is not quadratic in the fields: 

S = [ d d xdt \i4>(x,t)[d t g{x,t) -TV 2 g{x,t) - V ■ {gVfi)]+Tg(x,t)[Vi4){x,t)} 2 ) (36) 
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We are now interested in the linear response of a generic local functional of the density 
g to the external field p, which is defined as: 

/8A[g\(x,t)' 



Ra{x, x'\ t, t') 



\ 5p(x',t>) 



fi=0 



Integrating (|36|) by parts and taking the functional derivative of the generating 
functional with respect to p yields the field-theoretic expectation value: 

R A (x,x';t,t') = (A[g](x,t) V • [{V i$) g]{x' , f)) (37) 

The particular case A[g] = g makes explicit the distinction mentioned at the end 
of Section 12.2.21 between the actual response function (p(x,t)V ■ [(Vi(j>)g](x',t')) 
and the pseudo-response functions given by the propagators (i^>(x',t')p(x,t)) and 
(p(x,t)i(j)(x',t')). We now show how a symmetry of the action implies the fluctuation- 
dissipation theorem involving the proper response function defined in ([3j 



4-2. Symmetry of the action 

The action (fl4j) is invariant, modulo total derivatives, under the time-reversal 
transformation T, implemented as: 

t 

p(x,t) 
d t p{x, t) - 
i(f>(x,t) 



T : <^ 



T = -t 

p(x,r) 
-d T p(x,r) 
—i(j)(x, t) — In ^ ' 



(3f 



PO 



In general, the invariance of the action under a field transformation is not sufficient 
to guarantee that such a transformation is a symmetry: in fact, a non-linear field 
redefinition might lead to a non-trivial Jacobian. We now argue that, although the 
transformation ( l38i) is non-linear, its Jacobian is just a simple constant. Consider a 
finite system with discretized space: a field configuration is then defined by a 2iV-vector 



w] Pi, 



,Pn) an d the Jacobian matrix J associated to (!3~8l) is: 



J, 



d(Tv) t 1 


(a 


B 






D 



with: 



pr 1 














\ 



A — — Ijvxiv j B — - 

\ p^ / 

Since J is block triangular, |det(J)| = | det(A) det(D)| 
symmetry of the field theory. 



C = 



NxN 



D = I 



NxN 



± 11 and fl38l) is indeed a 
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The study of correlation functions is easier in Fourier space. The action (I14p can then 
be rewritten as: 



S 



k.t 



1 5F 
io{k.m n {-h:t) i$(k,t) / -Q(k,k\t)-—^- 



i0(fc,t) / Q(k,k',t)i$(k',t) 
Jk' 

where: J k = J j^y and J t = J dt. The integral kernel Q is defined as: 

Q(k,k',t) = Tk- k'p(-k- k',t) 
and Q~ l is such that: 

/ Q{k, k', t)Q" 1 (fc / , k", t) = 5{k - k") 
Jk' 

The time-reversal symmetry can now be implemented as: 



(39) 



T : I 



t — ► T — — t 

p(k,t) -> p(k,r) 
d t p(k,t) -> -d T p(k,r) 

4(k,t) -> -4(k,r)-J kll Q- 1 (k,k",t)d T p(-k",r) 
where the last line can also be expressed as: 

/ Q(k,k',t)4(k',t) -> - I Q(k,k',t)i(f>(k',T)-d T p(-k,T) 

Jk' Jk' 

Notice that f k , Q(k, k' ,t)i<fi(k f ,t) is nothing but the Fourier-transformed version of the 
composite operator V- [(Vi(f>)g] appearing in ( |371) . The FDT associated to any observable 
A is then: 

r 



(A(k,t)Q(k,k',t')i<j)(k',lf)) 



(40) 



(A(k, r)Q(k, k', T')4(k', r')} - d T , (A(k, r)p(-k, r')) 



In particular, we can consider A = p and t' < t so that one of the response terms 
vanishes because of causality. Assuming time-translation and time-reflection invariance, 
properties which are certainly true at equilibrium, we finally obtain the FDT: 



Tk ■ k'(p(k,t)p(-k - k',t')ij>(k',t')) = d t/ (p(k,t)p(-k,t')) 



(41) 



A few general properties useful in the following computations are: 

- All the FDT's written above hold only once the sum of connected and disconnected 
correlators is considered. 

- For the gas, any correlator containing more than one i<j> field vanishes. 

- (id>) EE 0. 
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- Each connected n-point function (n> 2) involving a field i(f) evaluated at the same 
time as another field p is identically zero. This can be easily argued from the 
integration domains in ( 1T6]) . (1TB]) . ( 125]) which vanish if m 2 ^ n = t\ whereas the 
integrand is a continuous function. 

Verifying ( )4T1) with A = p using the Feynman rules for the propagators is a simple 
exercise, yielding: 

- TPpoe-^^eit - t') = +Tk 2 p e- Tk ^ t '-^9(t' -t)- Tk 2 ^e-^^'lsgn^ - £)(42) 



When computing explicitly (14*0]) in the case A = pp or higher powers of p interesting 
cancellations occur between different disconnected correlators. Let us consider the 
simplest case (assuming t\ < t 2 < t 3 ): 

(p(k 3 , t 3 )p(k 2 , t 2 )Q(k 2 + k 3 , k', ti)i0(fc', £ x )) = <9 n (p(/c 3 , r 3 )p(k 2 , r 2 )p(-k 3 - k 4 , n)) 



After substituting (1391) and expanding in the non-vanishing disconnected parts, the 
formula above can be written in short-hand notation as: 



(43) 



T{k 2 + k 3 ) ■ k' ' \{p 3 p 2 i(f>i)c(pi) + W0i)c(P2Pi) c + 

(p2i4>l)c(p3pl}c + (P3#l)c(p2)(pi) + (P2i}l)e{p&){pi) 
= 9 n ((p 3 p 2 pi) c + (p 3 pi) c (p2) + (P2Pi)c(P3) 

In particular we now show that the following equality holds: 

T(k 2 + k 3 )-k'((p 3 p 2 i4> 1 } c (p 1 } + (p 3 i(j)i}c(p2Pi}c+ (p2i4>i)c(p3Pi}c) = d Tx (p 3 p 2 px) c (44) 



The remaining terms appearing in (1431) and not in (1441) contain 1 and 2-point functions 
only and are easily proven to be equal using (T42l) . A key observation is that conservation 
of momentum implies different values for k! in the l.h.s. of ( j4*4l) : k! = —k 2 — k 3 , k' = —k 3 
and k' = —k 2 respectively. Then, graphically: 

1 

,v2 A l 

Tp 



(k 2 + k 3 ) ■ (k 2 + k 3 ) ( 



2 3 



(fe + fc 3 ) -^3/ \ ~(k 2 + k 3 )-k 22 



2 3 

a 



(45) 



which is clearly true after simplifying the wedge graphs and considering the graphical 
rules of section 12.3.11 There are two points worth noting: first of all, the non-linear 
symmetry leads to cross-cancellations between disconnected correlation functions of 
different rank. Second, the particular dot-product structure of the vertex, together 
with conservation of momentum, once again produce interesting combinatorial factors 
which, in the particular case of (T4"5"l) . allow a single wedge graph to cancel two wedge 
graphs. A generalization of this kind of computations is not a hard exercise. 
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5. Conclusions and outlook 

As a complementary approach to the direct development of approximations for an inter- 
acting fluid, in this paper we have analyzed a non-interacting gas of Brownian particles; 
the system has been fully characterized with the computation of all the n-point func- 
tions. The exact correspondence between the particle and the field-theoretic formalism 
has been established by the non-perturbative calculation of the ground state of the field 
theory and the determination of the Poissonian probability distribution for the density 
field. A non-perturbative treatment of the gas- vertex is thus necessary not only to pre- 
serve the fluctuation-dissipation theorem in bare perturbation theory, as already argued 
in previous works [I], but also, and more fundamentally, to define the unperturbed 
ground state as the starting point for a perturbative expansion of physical observables 
in the inter-particle potential. 

At the fundamental level, we have shown that the system can be described exactly 
by a field theory and that no coarse-graining procedure is required to define the density 
field appearing in the path integral formulation. Hence, when considering the interact- 
ing theory, any divergence arising from the computation of loop diagrams cannot be 
ascribed to the effective nature of the field-theoretic description: instead it should be 
considered as a genuine feature of the problem and properly renormalized. Along the 
same lines, the particle origin of the field theory, expressed through the vertex, might 
determine and naturally constrain the renormalizability properties of the interacting 
theory. 

Finally, the pattern through which the FDT is verified on the correlation functions of 
the gas could shed light on its effect on the combinatorics of the interacting case, where 
the interest is in the study of high order correlators to define a diverging dynamical 
length scale [T2] . 
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Appendix A. Definitions and conventions 

Fourier Transforms For a finite system, in standard notation and using 5k as a 
shorthand notation for 8k,o- 

/*(*) = h f d d x e-* k - x f(x, t) ; /(*, t) = £ e^f k (t) 
^ v k 
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yy* ,x = W(aO ; / d d xe~ ikx = V8 k 

In the field-theoretic formalism there is no need to consider a finite number of particles 
and a finite volume as an intermediate step, the density being the only meaningful 
quantity. The formulae above are then modified as follows: 

/(*, t) = j d d x e- ik -*f(x, t) ; f(x, t) = j e *»/(fc, t) 

J ^~/ k ' X = ; / ddx = (?«) d 8{k) 

where all the integrations extends over the entire <i-dimensional space. 

Note that due to the different definitions for the Fourier transforms, the dimensions of 
the density field defined in (jSj) are different in Fourier space between the particle and the 
field-theoretic formalism. For particle systems: Qk=o{t) = h f v d d x g(x,t) = po with di- 
mensions of volume -1 ; in the field-theoretic formulation: g(k = 0, t) = j d x g(x, t) = N 
which is dimensionless. 



Averages and a useful identity 

- (A) v indicates the average of the noise-dependent observable A over the realizations 
of the noises rji, ... , r] N . 

- (A)j_p 4 indicates instead the average over the initial positions of the particles. It is 
defined in such a way that (l)/.p. = 1: 



(A) 



LP. 



I J d^O)...^ J d d x N (0) A 



In the particle formalism, with the above conventions for the Fourier transforms, identity 
(tO|) follows straightforward result: 



fci+...+fc a _i+fc a+ i+...+fe n "fe a 



if all j's are the same 
if all fs but j a are the same 

(A.l) 



hi-Aa-iha+x-hb-Ab+i-Vka+h if all fa differ but 

Ja 3b 



k 4l • • • h„ 



if all j's differ 



Appendix B. Combinatorial coefficients 

Appendix B.l. Factorials 



Theorem: the coefficient TZi (defined in section 2.3.4}) associated to a vertex with I 
dangling lines is TZi = (I — 1)!. 



We have verified in the text that the relation is true for some I: TL2 = 1 = (2 — 1)!, as 
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shown in (I16p and in the comment below it. Formula (1191) implies 71% = 2 = (3 
From section [2.3.31 we also know that TZ^ = 6 = (4 — 1)!. 



Assuming that the relation is valid for vertices having up to / dangling, we now prove it 
for a vertex with I + 1 dangling lines. Let us consider one such a vertex, with the upper 
edge labeled by I + 1. 



K, 




+ 





(B.l) 



+ 



+ 




In a Dyson-Schwinger fashion the earliest vertex is pulled out of the blob. Conversely, 
when contracting the edges generated by the extraction of the vertex, the graphs on the 
r.h.s. generate the graph on the l.h.s.. Note that in (IB.lj) any vertex with two dangling 
lines must be a bare vertex since no loop-graphs contribute to the physical amplitudes. 
All binary (non-planar) trees are generated by a recursive application of (IB.lj) . They all 
appear with coefficient 1 except for one particular case: when I is even, the graph on the 
r.h.s. for which r = s = 1/2 appears only once and the factor of 1/2 out of the braces in 
(IB.lj) is not simplified. However such a graph has a left-right symmetry and because of 
the very same symmetry a factor 2 will arise in formula ( IB. 31) since the partition with 
r elements and the partition with s elements are equivalent. This cancellation will be 
tacitly assumed in the argument following ( IB. 31) . 



Every term in the r.h.s of ( IB. II) represents a partition of the / dangling lines into two 
groups with r and s dangling lines each, with the constraint I = r + s. For a generic 
partition, the Feynman rule associated to the extracted vertex reads: 

- 2 ( fc i + • • • + K) ■ (k r +i + ... + k) (B.2) 



TI 



The X/r/ * s over ^ ne C r ) possible partitions depicted in (IB. II) . Applying conservation of 
momentum, we can rewrite ( IB. 21) as: 

2 J2( kl + ■■■ + k r )- {ki+i + h + . . . + k T ) = 



(B.3) 



Tl 



TI 



(kl + ... + k 2 r ) + (J2 k i- k i) + k i+i -{ki + ... + k r 

The sum over all TI labelings can be easily computed by the following argument: in the 
sum each kf appears (' jl) times, since once ki is fixed to be in the group of r elements, 
the other r — 1 elements can be chosen among / — 1 labels. Similarly, each ordered pair 
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kikj in the inner sum appears ( _n) times, since once ki and kj are fixed to be in the 
group of r elements, the other r — 2 elements can be chosen among I — 2 labels. 

= (lz i) (*? + ... + *?) + (' : i) e ** • + (' : i) k W . & + ... + *o ( B . 4 ) 

= f : i) <*? + . . . + *?> + g : 2 2 ) e *. • (-*, - - (' : i) ^ 
= (' : i) (-^ 2 +1 + ^ + . . . + k f ) + g - 2 2 ) E(-^ 2 + fc m) 



/-2 
r — 1 



"^f+i + ^i + • • • + k i 



where in the second and third steps, conservation of momentum has been used. 
Summing over all possible partitions of labels into two blobs and associating to each 
blob with n dangling lines the appropriate lZ n — (n — 1)!, we finally obtain: 

Ef^V-W-r-l)! = (i-2)!^l = (1-1)! (B.5) 

r=l ^ ' r=l 

Appendix B.2. Generalized Stirling numbers 

Theorem: the coefficient C n (defined in section ^. 4 -4\ ) associated to the c-vertex with n 
dangling lines appearing in the computation of c-functions is: C n = S n -2, where S n is 
the sequence of generalized Stirling numbers (Sloane A001705, /77]/ ). 

We can write a Dyson- Schwinger equation similar to (IB. II) . with the important dif- 
ference that the c-vertex needs to have 3 dangling blobs: the first blob with 1 dangling 
lines, the second one with m, and the third one with n — I — m. The factor 1/2! in 



( IB. II) becomes 1/3! taking into account the possible permutations of blobs. The rig- 
orous argument about the cancellation of symmetry factors for symmetric diagrams is 
analogous to the one which followed equation (IB. lh . Formula (1B.2I) is replaced by the 
Feynman rule for the c-vertex obtained in (126]) : 

T\ E [( kl + - + ^ + ( kl+1 + - + kl+m ^ + + - + fc «) 2 ( R6 ) 

' TI 

In the sum over TI relabelings each kf appears once for each of the possible (") ( n ~ l ) 
labelings, since ki has to be in one of the three groups in (1B.6j) . If we ask an ordered 
pair k^kj to appear in the blob with 1 dangling lines, this happens (™_|) (^ n 2 \^ 1 2 ^) 
times, since once ki and kj are fixed to be in the group of 1 elements, the other 1 — 2 
labels can be chosen among n — 2; the other labels are then partitioned into blobs with 
m and n — 1 — m dangling lines respectively. Focusing instead on the blobs with m 
and n — I — m a given ordered pair appears [ n ~V\ ( n ~i m ) an d (' n_2 )~^~ 2 -') ^ mes 
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respectively. Finally, to each blob with r dangling lines we associate a factor (r — 1)! as 
derived in the previous section. The sum over TI labelings can now be written as: 

1 n—2 n—l—1 

^jE E(*-l)!(m-l)!(n-i-m-l)!. 

1=1 m =l 

'n\ (n — t 
I / \ m 



E*. 2 



i=l 



n — 2\ fn — I 



I - 2 



m 



+ 



n 



m 



n — m 
I 



+ 



n-2 



n — I — m — 2 



I + m N 
I 



(B.7) 



Exploiting once more conservation of momentum: Yl^j ' kj = ~Yli=i^i- ^ e can 
then factor out Y17=i which always cancels the denominator of the rational function 
generated by the c-vertex and will be dropped in the following. Simplifying the binomial 
coefficients we obtain: 

n-2 n-l-1 



sEE <»-*>'( 



n — I 



+ 



n — m 



m{n — I — m) l(n — I — m) 



+ 



I + m 
Im 



1=1 m =l 

Moreover, the three terms in parentheses clearly give the same result when summed 
over I and m: for simplicity we focus only on the third one and multiply by 3. We thus 
have, exchanging the order of the summations in the second term: 

n—2 n—l — 1 , n—2 n—m—1 



>-KEE^EEy) 

1=1 m=l m=l 1=1 

n-2 n-2 

-{n-2)\(^Hn-i-i + Y, H ^-i) = S n 



1=1 m=l 

where by H n we have indicated the harmonic numbers. 



Appendix C. Graphology 

Appendix C.l. Correlation functions are represented by ladder graphs 
Theorem: the analytic expression associated to an n-point ladder graph, namely: 



exp 



kl(t 2 -h) 



(ki + ... + k n _i) 2 (t n - t n _i) 



(CI) 



is the same as the analytic formula (TJ|) for the n-point connected correlation function 
computed in the particle formalism: 



n 

{g(k 1 ,t 1 )...Q(k n ,t n )) ri = expfc^ki-kjfa-tij) 

i<j 



(C.2) 



Examples of ladder graphs are drawn in figure (133]) . and the absolute value of the time 
differences in flC.2j) has been removed considering the convention ti < . . . < t n . 
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Let us define the variables k s as kj = Y2i=i k%- Hence we can rewrite fIC.lj) as: 



n-1 



exp ( ~ Yl ^(Wi - tm)) (C.3) 



m=l 



We will now prove the theorem, showing that the exponent in ( 1C.2I) can be rewritten as 

the exponent in (IC.3I) . 

Expanding (tj — ti) in a telescopic sum: 

j-i 

tj — ti = {tj — tj-i) + (tj-i — tj-2) + • • • + (t i+ i — ti) = ^J(t m+ i — t rn ) 

m=i 

we can write: 

n j-1 n j-1 j'-l 

/] fa ■ kj{ti — tj) = ki ■ kj{t m — t m+ i) 

j'=l i=l j=l i=l m=i 

Exchanging the order of the summations we obtain: 

n j—l j—l n—1 m n 

EEE — EE E (C4) 

i=l i=l m=i m=l i=l i=m+l 

and looking at the innermost summation, conservation of momentum implies: 

n m 

(y ^ j kj^ ■ ki(t m — t m+ i) = — ^ ^ kj ■ ki(t m — t m+ i) (C-5) 

j=m+l j=l 

Applying (1C.4j) . (1C.5[) and the definition of k it is now easy to pass from the exponent 
in ( 1C.2I) to the exponent in ( 1C.3I) : 

n j — l j — l n—1 m n 

h • kj(t m - t m+ i) =yy y h ■ kj(t m - t m+1 ) (C.6) 

i=l i=l m=i m=l i=l j=m+l 

n—1 m m n—1 

= ~ ^ ^ ^ ^ kj ■ ki(t m — t m+ i) = — ^ k^itm — t m+ i) (C-7) 

m=l i=l j=l m=l 

Appendix C.2. Cancellations among r-graphs with a stalk 

Theorem: in the sum over topologically inequivalent (TI) labelings, r-graphs with a 
stalk cancel one another, leaving only one n-ladder graph. 

Examples of r-graphs with a stalk appear in the first line of (!34j) . Let us assume that 
one of the r-vertices of a graph has n > 2 dangling lines. The combinatorial factor 
associated to such vertex is TZ n — (n — 1)!. In the sum over TI labelings, graphs with 
such a vertex can be obtained in several ways: either from a labeling of the graph itself 
or from a labeling of a "less contracted" graph which induce one or more contractions 
generating the vertex. Since n > 2 we know that at least one contraction has occurred. 
We can now work backwards and undo the contraction to obtain all the possible graphs 
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which could have generated the original graph. From the group of n labels we can pull 
out a group of k labels, leaving a vertex with n — k + 1 dangling lines in such a way 
that the largest label is in the extracted group. The rules in section 12.3.41 imply then 
that the edge we have just created must indeed be contracted. The combinatorial factor 
associated to the two vertices just created are {k — 1)! and {n — k)\ respectively. Note 
that since the largest label has to be in the group with k elements, among the n — 1 
labels left we have to choose k — 1 of them. Summing over all possibilities: 



k=2 v 7 

Remember now that each contraction brings a factor of (—1) so that all the graphs we 
generated with one contraction less than X appear with a relative minus sign. But 
(1C.8I) implies that the numerical value of the coefficient is the same and the sum over 
the whole group of graphs vanishes. 

On the other side ladder diagrams survive. In fact, in order to create a ladder, from a 
blob with n labels we pull out recursively a group of n — 1 labels which do not contain 
the earliest time among the n labels. Then, the argument above leading to the graph 
cancellations does not apply. 

Appendix C.3. c-graphs cancel r-graphs without a stalk 

Theorem: in the sum over topologically inequivalent (TI) labelings, c-graphs completely 
cancel r-graphs without a stalk. 

Examples of r-graphs without a stalk and c-graphs appear in the second and third 
line of ( 134"|) respectively. 

Let us consider an n + 1 point function, r-graphs without a stalk have n leaves and c- 
graphs have n + 1 leaves. It is useful to note that for r-graphs the label '1' is constrained 
by causality to appear on the upper vertex (root), as a consequence of the contraction 
of the stalk. For c-graphs, instead, the c- vertex itself is always labeled by '1' ('1' being 
the smallest among all the descendant labels), another label '1' has to be placed on 
some leaf, and because of this some edges need to be contracted according to the rules 
discussed in section 12. 3.1[ The key idea of the proof is that for an r-graph the number of 
edges dangling from the root is independent of the labeling, whereas for c-graphs edge 
contractions cause this number to vary. 

Let us consider an r-graph X whose root has / dangling lines with the associated combi- 
natorial factor (I — 1)!. It is possible that particular labelings of some c-graphs Y which 
initially have 3, . . . , / + 1 dangling lines produce X after operating the contractions 
induced by the labelings. Those graphs will be multiplied by the generalized Stirling 
numbers Si, ... , Sj_i. In particular, we can construct all the graphs Y starting from 
a c-vertex with / + 1 dangling lines and "pulling down" groups of edges with the con- 
straints that a c-vertex must have at least three dangling lines and that '1' labels one 
of the lowest leaves in such a way to induce a contraction leading back to a root with I 




(C.8) 
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dangling lines. 

An example should make things clear. One of the r-graphs without a stalk contributing 
to the computation of the 5-point function, with the appropriate factor of (4 — 1)! = 6, 
is: 



X = Q //\\ (C.9) 

Among the c-graphs, the ones which follow generate X if, as explained above, '1' labels 
one of the lowest leaves. The coefficient in front of each graph is made up of: a sign 
which keeps track of the number of contracted edges; the generalized Stirling number 
associated to the c-vertex; the product of coefficients associated to the underlying r- 
vertices; the number of labelings which leave '1' to label one of the lowest leaves. 



Y's = -26-1-1 Jb\ + 5 • 1 • r\ 



or equivalently: 

-S 3 + S 2 Q + ft Q) 2! + ft (f\ (fj = -(4 - 1)! (CIO) 

We see that (I — 1)! is exactly canceled by a linear combination of the generalized Stirling 
numbers ft, ... , ft_i. In this case / = 4. 

Working out explicitly more complicated graphs through the "pulling down" proce- 
dure, it is not hard to write the general form for the l.h.s. of ( IC.lOp . which is, modulo 
a sign: 

1-2 / , x 1-2 k 2 -l 



*-> - E>—(*> ! ^mjKBft-w 
EEE^(0(SO».-*- + - 

I i ,-. Q l-^. O In, 1 \ / \ / \ / 



k 3 =3 fc 2 =2fci=l 

1-2 ^ 1-2 ^ k 2 -l 

Si-i - E^-^(rr^)T + E^-*-* jrrj^i E 1 

fcl=l V 1J k 2 =2 V ZJ fci=l 

1-2 ,j fca-lfa-l 



^ + ee<-^.-^C':;) 



j=i kj=j 

1-2 ,, k 



5,-1 - = (n-1)! 
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